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Abstract 

Background: Strawberries {Fragaria ananasso) reproduce asexually through stolons, which have strong tendencies 
to form adventitious roots at their second node. Understanding how the development of the proximal (1-1) and 
distal (I-2) internodes of stolons differ should facilitate nursery cultivation of strawberries. 

Results: Herein, we compared the proteomic profiles of the strawberry stolon 1-1 and I-2 internodes. Proteins 
extracted from the internodes were separated by two-dimensional gel electrophoresis, and 164 1-1 protein spots 
and 200 I-2 protein spots were examined further. Using mass spectrometry and database searches, 38 1-1 and 52 
I-2 proteins were identified and categorized (8 and 10 groups, respectively) according to their cellular 
compartmentalization and functionality. Many of the identified proteins are enzymes necessary for carbohydrate 
metabolism and photosynthesis. Furthermore, identification of proteins that interact revealed that many of the I-2 
proteins form a dynamic network during development. Finally, given our results, we present a mechanistic scheme 
for adventitious root formation of new clonal plants at the second node. 

Conclusions: Comparative proteomic analysis of 1-1 and I-2 proteins revealed that the ubiquitin-proteasome 
pathway and sugar-hormone pathways might be important during adventitious root formation at the second node 
of new clonal plants. 



Background 

Several stoloniferous species produce long, sympodial 
stolons with rooted rosettes (ramets) at their nodes 
[1,2]. For the garden strawberry {Fragaria ananassa), 
the mother plant forms plantlets on stolons during 
spring growth (Figure 1). The first stolon originates 
from an auxiliary leaf bud produced in the central 
crown and commonly contains only two nodes. The 
regions along the stolon and between the plant and the 
first node and the first and second nodes are the 1-1 
(proximal) and 1-2 (distal) internodes, respectively. 
Although stolon growth requires internode elongation, 
the fates of the two nodes are dissimilar. First, 1-1 elon- 
gates and terminates at the first node, which is nonpro- 
ductive, then 1-2 elongates and terminates at the second 
node, which forms the main crown of the clonal plant. 
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Many studies have assessed the relationships between 
plant genotype and phenotype [3] using morphological 
differences caused by loss of function or altered expres- 
sion of a single gene. To fully understand the function 
of a gene, however, the expressed protein must be char- 
acterized. Proteomics investigates the synthesis, turn- 
over, and modification of proteins so that gene function 
and genotypes can be understood [4]. For example, a 
proteomic study of ripening strawberry fruit from plants 
of different genotypes identified constitutively and dif- 
ferentially expressed proteins that probably control the 
quality of the fruit [5]. Proteomics may, therefore, be 
used to address biochemical and physiological aspects of 
plant morphologies. Such approaches are increasingly 
used to elucidate the biochemistry and physiology of 
model species [6,7]. Yet proteomics is limited by the 
ability to identify proteins, which relies on the availabil- 
ity of sequence data. For plant species with unsequenced 
genomes, proteomics can still be applied but the num- 
ber of identified proteins is usually smaller because their 
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Figure 1 The strawberry stolon. As shown, a stolon shoots away 
from the base of a strawberry plant. A clone is formed at a variable 
distance away from the parent at the second node concomitant 
with adventitious root formation. 



identification relies on homology with proteins of other 
species. For example, although F. ananassa sequence 
data is very limited, two proteomic techniques, namely 
two-dimensional gel electrophoresis (2-DE) and mass 
spectrometry (MS), have been used to identify an 
F. ananassa protein homologous to the birch pollen 
allergen Bet v 1 [8]. 

Elongation of the strawberry stolon is considered to be 
the result of cell division and cell expansion [9], but lit- 
tle is known about how the 1-1 and 1-2 internodes 
develop. Such knowledge is required, however, if we are 
to improve the cultivation of nursery strawberries and 
understand in greater depth how clonal multiplication 
occurs. 

For the work reported herein, we compared the pro- 
teomes of the F. ananassa 1-1 and 1-2 internodes to elu- 
cidate the differences in their growth and functional 
characteristics and establish reference maps by identify- 
ing the protein spots of their 2-DE maps in conjunction 
with MS peptide mass determination and database 
searches. We identified isoforms of several proteins and 
present a detailed analysis of the two proteomes, which 
allows us to begin to explore the different developmen- 
tal mechanisms of the 1-1 and 1-2 internodes. The refer- 
ence maps should be useful for investigation of 
strawberry physiology and for monitoring changes in 
protein expression in strawberry stolons in response to 
biotic and abiotic stresses. 

Materials and methods 

Plant material 

The F. ananassa cultivar Hongjia was obtained from a 
nursery at the Hangzhou Academy of Agricultural 
Sciences, Zhejiang, China. The plants were grown in a 
tunnel greenhouse with a 10-h light/ 14-h dark cycle, a 
30°C-day/26°C-night temperature cycle, 150 umol m" 2 s" 1 
light intensity, and a relative humidity of 60%. Plants 



were watered regularly and provided adequate nutrients. 
In the morning, after transplantation, stolons that had 
formed from the leaf buds of three-month-old plants 
were collected and cut to isolate the 1-1 and 1-2 inter- 
nodes (Figure 1). Daughter plants that had formed at the 
apices of stolons had not been allowed to root. For each 
experiment, 50 stolons were randomly chosen and 
removed from four or five plants. The 1-1 and the 1-2 
internodes were pooled separately, rinsed with water to 
remove contaminants, quickly dried with paper towels, 
frozen in liquid nitrogen, and stored at -80°C prior to 
protein extraction. 

Protein extraction 

Proteins were extracted using acetone and trichloroacetic 
acid method. A portion (2 g) of each internode sample 
was pulverized with a pestle in a mortar that contained 
liquid nitrogen and then homogenized in 10 mL of 10% 
(w/v) trichloroacetic acid, 0.07% (v/v) 2-sulfanylethanol 
in acetone. Total protein was precipitated for 1 h or over- 
night at -20°C. The extracts were each centrifuged at 
13000 x g for 20 min at 4°C. The pellets were washed 
three times with 0.07% (v/v) 2-sulfanylethanol in acetone, 
and vacuum dried for 30 min. The dried powders 
(30 mg) were each resuspended in 500 uL of 7 M urea, 
2 M thiourea, 4% 3-[(3-cholamidopropyl) dimethylam- 
monio]-l-propanesulfonate (CHAPS), 0.75% dithiothrei- 
tol (DTT), 0.5% Biolyte (pH 3.0-10.0, Bio-Rad), 1 mM 
phenylmethanesulfonyl fluoride, and then shaken vigor- 
ously for 1 h at room temperature. Insoluble material 
was removed by centrifugation at 13000 x g for 15 min at 
20°C. At least three replicates were prepared. Protein 
concentrations were determined using Bio-Rad Protein 
Assay kit reagents (standard Bradford method) with 
bovine serum albumin as the calibration standard [10]. 

Two-dimensional gel electrophoresis 

Each sample contained 300 ug protein in 350 uL of 8 M 
urea, 2 M thiourea, 2% CHAPS, 0.5% Biolyte (pH 3-10), 
0.75% M DTT, 0.002% Bromophenol Blue. Each sample 
was each loaded onto a 17-cm immobilized pH (3-10) 
gradient strip (Bio-Rad). The strips were rehydrated for 
12 h at 50 V. Isoelectric focusing used a linear ramp 
from 0 to 250 V in 15 min, a linear ramp from 250 to 
10000 V in 1 h, and 10000 V for 5 h, all at 20°C. After 
isoelectric focusing, the strips were equilibrated in 
50 mM Tris-HCl, pH 8.8, 6 M urea, 20% glycerol, 2% 
sodium dodecyl sulfate (SDS), 2% DTT, and then in a 
solution of the same composition that also contained 
2.5% (w/v) iodoacetamide, (the time of each incubation 
was 15 min). The strips were then each placed onto a 
1-mm-thick SDS (12.5% (w/v)) polyacrylamide gel and 
sealed with 1% (w/v) agarose. Electrophoresis was 
carried out in a Bio-Rad PROTEAN apparatus at 
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24 mA/geL The gels were stained using a modified sil- 
ver-staining method that is compatible with MS [11]. 
Image analysis was subsequently performed. These pro- 
cedures were replicated three times. 

Image acquisition and analysis 

The three replicates of the 1-1 and 1-2 2-DE gels were 
scanned using a calibrated densitometer (GS-800, Bio- 
Rad), and the spot patterns were characterized using 
PDQuest software (ver. 8.0.1, Bio-Rad). Image analysis 
steps included image filtration, spot detection and mea- 
surement, background subtraction, and spot matching. 
One 1-1 gel served as the reference, and the spots of the 
other five gels were referenced to it. Initially, spots were 
automatically matched, and the positions of unmatched 
spots were then manually determined. The molecular 
mass (kDa) of each protein was estimated by compari- 
son with those of a standard marker set, and the isoelec- 
tric points (pis) were determined by the spot positions 
along the immobilized pH gradient strips. 

In-gel protein digestion and mass spectrometry 

The silver-stained protein spots were manually excised 
from the gels, and each was placed into a well of a 
96-well microplate. The gel pieces were destained in a 
solution prepared from a 1:1 (v/v) mixture of 30 mM 
potassium ferricyanide and 100 mM sodium thiosulfate 
at room temperature for 10 min, vortexed until 
destained, washed three times with 300 [iL of Milli-Q 
water (each time for 5 min) and dehydrated in 150 [iL 
of acetonitrile. Then the gel samples were swollen in 
50 mM NH4HCO3 containing 12.5 ng/(iL trypsin 
(Sigma, Cat. No. 089K6048) at 4°C for 30 min, and at 
37°C for longer than 12 h. For each digest, the peptides 
were extracted from the gels twice with 5% trifluoroace- 
tic/50% acetonitrile at room temperature, resuspended 
in 0.7 \iL of 0.2 M a-cyano-4-hydroxy-cinnamic acid 
(Sigma) in 0.1% trifluoroacetic/50% acetonitrile, and 
allowed to dry under a stream of nitrogen. The 
extracted peptides were subjected to matrix-assisted 
laser desorption/ionization time -of- flight MS (4800 Pro- 
teomic Analyzer Applied Biosystems). Proteins were 
identified using the Peptide Mass Fingerprinting module 
of Mascot (Matrix Science) and the experimental masses 
of the peptides. We searched the Swissprot database in 
September 2010 (version 20100906, which included 
519348 sequences and 183273162 residues) for proteins 
from Viridiplantae (green plants, 29439 sequences). One 
missed cleavage per peptide was allowed, and a mass 
tolerance of 50-150 ppm was used. Carbamidomethyla- 
tion of cysteine was set as a fixed modification, and oxi- 
dation of methionine was allowed. Identified proteins 
with a peptide mass fingerprint were denoted as having 
an unambiguous identification by the following criteria: 



(1) at least five different predicted peptide masses were 
needed to match the observed masses for an identifica- 
tion to be considered valid; and (2) protein scores needed 
to have >57 identity for Swissprot database (p < 0.05). 

Gene ontology (GO) annotation and protein classification 

The UniProt database (http://www.uniprot.org) was 
searched to determine the functions of the identified pro- 
teins. Three independent ontological sets in the Viridi- 
plantae taxonomic database were used to annotate and 
group the proteins according to biological processes, 
molecular function, and cellular compartmentalization. 

Protein-protein interactions (PPIs) 

PPIs were predicted by Cytoprophet, which is a plug-in 
of Cytoscape. We sent all the data for the identified 
proteins in 1-1 and 1-2 to the Cytoprophet server along 
with the set cover approach maximum specificity set 
cover (MSSC). The algorithm was used to predict the 
interaction network(s). 

Results and discussion 

Proteome analysis and protein identification 

We investigated the differences between the 1-1 and 1-2 
protein profiles. More than 503 1-1 and 1127 1-2 protein 
spots were reproducibly detected (Figure 2). The pis of 
the protein spots ranged from 3.5 to 9.3, and the mole- 
cular masses ranged from 7.1 to 60.2 kDa. 

For protein identification, the peptide mass finger- 
printing data were used in conjunction with a search of 
the complete Swissprot Viridiplantae taxonomic data- 
base as only 1017 F. ananassa protein sequences were 
available therein. Respectively, 164 and 200 spots from 
the 1-1 and 1-2 gels were selected for MS, and 38 and 
52 proteins were identified, i.e., -25% of the total in 
each case (Tables 1 and 2). No unambiguous matches 
were made for the other 274 proteins, probably because 
the proteins were not included in the database or 
because a protein spot contained more than one protein. 
For the 90 identified proteins, 5 (6.5%) were common to 
both proteomes, which indicated that the two internodes 
had some proteins (probably housekeeping proteins) in 
common. Certain proteins (35% of the gene products) 
were specific to 1-1, whereas 58.5% were specific to 1-2 
(Figure 3). More than 50% of the proteins could be cor- 
related with annotated proteins from at least one dicoty- 
ledon species (Figure 4). Of the F. ananassa proteins 
identified in our study, 13 1-1 and 25 1-2 proteins were 
matched to Arabidopsis thaliana proteins, and others 
were matched to Oryza sativa proteins (two from 1-1, six 
from 1-2) and Zea mays proteins (three from 1-1, two 
from 1-2). Only five 1-1 proteins and three 1-2 proteins 
could be matched to those found in the F, ananassa data- 
base, possibly because the number of annotated proteins 
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Figure 2 Two-dimensional SDS-PAGE gels of the 1-1 and 1-2 proteomes. Proteins (300 pig) in 1-1 and I-2 extracts were separated, in the first 

dimension by isoelectric focusing (pH 3-10) and in the second dimension by SDS-PAGE through 12.5% acrylamide gels. Proteins were visualized 

by silver staining. Circled proteins were identified by matrix-assisted laser desorption/ionization time-of-flight MS and database searches. 
I J 



contained in the database is relatively small compared 
with that for the A. thaliana database. 

Protein isoforms and subunits of protein complexes 

Isoforms for malate dehydrogenase (spots 15, 90, and 
91), glucan endo-l,3-P-glucosidase 6 (spots 29 and 31), 
triosephosphate isomerase (spots 69 and 72), oxygen- 
evolving enhancer protein 1 (spots 126 and 127), and a 
putative F-box/kelch-repeat protein (spots 62 and 151) 
were found in the 1-1 proteome, and isoforms for malate 
dehydrogenase (spots 129 and 141), U-box domain-con- 
taining protein 34 (spots 153 and 175), and the probable 
protein ABIL5 (spots 12 and 48) were found for 1-2 
(Tables 1 and 2). Therefore, respectively, -29% and 
-12% of the identified 1-1 and 1-2 proteins exist as iso- 
forms. Previous studies reported that -70% of maize or 
Arabidopsis proteins exist as isoforms [12,13]. The 
experimental masses and/or pis of certain proteins 
differed from their theoretical values, possibly owing to 
co-translational and/or post-translational modification 
(e.g., glycosylation, phosphorylation, and/or proteolysis), 
translation from alternatively spliced mRNAs [14-17], or 
post-translational modification by a nonprotein compo- 
nents) [18]. 

Of the 90 identified proteins, 16 (more than 17%) 
form complexes [19]. For example, two 1-1 spots corre- 
sponded to the a and P subunits of chloroplastic ATP 
synthase. The P-subunit of the DNA-directed RNA 
polymerase subunit was detected in both the 1-1 and 1-2 
proteomes. For some complexes, all component subunits 
were detected; for instance, both the small and large 



subunits of ribulose-bisphosphate carboxylase/oxygenase 
(RuBisCo) were identified in the 1-1 proteome. The fact 
that several complexes were identified will facilitate 
future studies of how complex formation is regulated. 
The components of a protein complex would be 
expected to be coordinately expressed and regulated to 
keep a system in functional balance. Thus, our 2-DE 
reference maps can be used to compare the levels of 
functionally related polypeptides (isoforms and subunits 
of complexes) and may provide insights into protein 
function and participation in molecular networks. 

Functional classification and subcellular localization of the 
identified proteins 

Identification of proteins that are differentially expressed 
in 1-1 and 1-2 is important for our understanding of 
internodal development and differentiation. The identi- 
fied proteins were grouped according to their biological 
processes, i.e., photosynthesis, protein synthesis, protein 
folding, transcription, transport, stress, and development, 
and cellular locations using the GO annotation in the 
Viridiplantae taxonomic databases (Figures 5 and 6). 
Proteins involved in metabolism and photosynthesis 
accounted for 32% and 18%, respectively, of the 1-1 pro- 
teins. The numbers of 1-1 proteins involved in transcrip- 
tion (13%) and energy production (8%) were also 
substantial. Interestingly, for 1-2, although still the largest 
two groups, metabolic proteins accounted for only 15% 
and those involved in photosynthesis for 16%. 1-2 pro- 
teins involved in protein synthesis and DNA transcrip- 
tion each accounted for 12% of the total. The function 
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Table 1 Identified proteins from the strawberry stolon 1-1 internode 

Spot Protein Accession Molecular Reference Theoretical Experimental Score f SC 9 Matched/ 

No.° Number 6 Function c Organism kDa/pl d kDa/pl e (%) Unmatched 

queries' 1 





Metabolism 


















15 


Malate dehydrogenase 


P83373 


L-malate 


Fragario 


35.8/8.7 


13.5/5.3 


99 


41 


9/43 








dehydrogenase 


ananasso 


















activity 














29 


Glucan endo-1,3-beta- 


Q93Z08 


glucan endo-1,3- 


Arobidopsis 


52.6/5.6 


20.3/5.8 


58 


27 


6/35 




glucosidase 6 




beta-D-glucosidase 


thaliana 


















activity 














31 


Glucan endo-1,3-beta- 


Q93Z08 


glucan endo-1,3- 


Arabidopsis 


52.6/5.6 


20.4/5.0 


63 


21 


5/42 




glucosidase 6 




beta-D-glucosidase 


thaliona 


















activity 














40 


Uricase 


004420 


urate oxidase 


Arobidopsis 


35.0/8.6 


23.4/6.1 


58 


38 


10/107 








activity 


thaliana 












42 


GDSL esterase/lipase 


023469 


hydrolase activity, 


Arabidopsis 


26.7/8.9 


23.2/7.8 


57 


37 


6/36 




At4g 16220 




acting on ester 


thaliana 


















bonds 














90 


Malate dehydrogenase, 


P83373 


L-malate 


Fragaria 


35.8/8.7 


40.3/8.2 


122 


51 


10/50 




mitochondrial 




dehydrogenase 


ananassa 


















activity 














91 


Malate dehydrogenase, 


P83373 


L-malate 


Fragaria 


35.8/8.7 


43.3/8.3 


100 


50 


9/45 




mitochondrial 




HphvHmnpna^p 


ananassa 


















activity 














1 1 1 


^-arlpnn^x/l mpthinni np 


ORW1Y4 


rnptDl inn hinrlinn 


Phnspnli k 
r / iujcuiuj 


43.5/5.6 


47.4/5.8 


74 


24 


7/9 




synthase 






lunatus 












118 


2,3-bisphosphoglycerate- 


P30792 


manganese ion 


Zea mays 


60.7/5.2 


41 .3/4.9 


59 


24 


8/27 




independent 




binding 
















phosphoglycerate mutase 


















133 


Fructokinase-1 


A2WXV8 


ATP binding 


Oryza sativa 


34878/5.1 


36.1/3.8 


76 


23 


6/15 


142 


S-adenosylmethionine 


A9PHC5 


metal ion binding 


Populus 


42999/5.7 


46.5/6.3 


74 


18 


5/2 




synthase 4 






trichocarpa 












1 62 


Phosphoglucomutase 


DOQ on/i 
r93oU4 


magnesium ion 


Zea mays 


D3ZOD/b.4 


60.2/6.1 


63 


1 7 


9/30 








binding 
















Energy 


















1 52 


ATP ^x/nth^KP ihi mit 


OQMRRQ 

v_/ J\ V 1 1 \\ \J 


ATP hinHinn 




53.8/5.2 


45.3/4.8 


76 


27 


8/1 9 




beta, chloroplastic 






schreberi 












1 53 


ATP c\/nthaco ci ihi mit hota 
nl r by 1 III Icoc bUUUIIIl Ufcrld 


nm rsq 


ATP hinrlinn 
r\ \ r Ul 1 IUII iy 


r~)r\/vn cnti\/n 
\JlyZ.U jULIVU 


590/5 9 


48 1/5 4 


1 34 


44 


1 9/45 


1 


ATP synthase subunit beta 


vjOvjDr Z 


a i r uiriGing 


Castanea sativa 


JJ.O/J.J 




1 HQ 
I Uo 


D / 


I J/ZD 


156 


ATP synthase subunit 


A4QJA0 


ATP binding 


Aethionema 


55.3/5.2 


58.5/3.5 


114 


22 


10/13 




alpha 






cordifolium 














Photosynthesis 


















2 


Ribulose bisphosphate 


P 10795 


monooxygenase 


Arabidopsis 


20.5/7.6 


10.7/7.8 


62 


43 


7/101 




carboxylase small chain 
1 A 




activity 


thaliana 












37 


Ribulose bisphosphate 


098681 


ribulose- 


Zamioculcas 


49.9/6.3 


21.4/4.9 


79 


25 


10/78 




carboxylase large chain 




bisphosphate 


zamiifolia 














(Fragment) 




carboxylase activity 














46 


Ribulose bisphosphate 


P28439 


magnesium ion 


Pelargonium 


53.3/6.3 


23.7/7.5 


62 


32 


13/104 




carboxylase large chain 




binding 


hortorum 












69 


Triosephosphate 


Q9M4S8 


triose-phosphate 


Fragaria 


33.7/7.6 


30.0/5.1 


152 


55 


19/112 




isomerase 




isomerase activity 


ananassa 












72 


Triosephosphate 


Q9M4S8 


triose-phosphate 


Fragaria 


33.7/7.6 


31.3/5.6 


126 


59 


16/77 




isomerase 




isomerase activity 


ananassa 












126 


Oxygen-evolving 


P26320 


calcium ion binding 


Solan urn 


35.6/5.8 


35.6/4.4 


111 


33 


8/8 




enhancer protein 1 






tuberosum 












127 


Oxygen-evolving 


P26320 


calcium ion binding 


Solan urn 


35.6/5.8 


34.9/4.8 


113 


28 


9/21 




enhancer protein 1 






tuberosum 
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Table 1 Identified proteins from the strawberry stolon 1-1 internode (Continued) 





Transcription 


















6 


Histone acetyltransferase 
GCN5 


Q9AR19 


protein binding 


Arobidopsis 
thaliana 


63.5/6.0 


14.2/4.1 


59 


29 


11/66 


18 


Pentatricopeptide repeat- 


Q9SF38 


nucleotide binding 


Arabidopsis 


84.4/6.5 


1 1 .2/6.2 


80 


23 


14/49 




containing protein 






thaliana 














At3g09650 


















27 


Transcription factor 
DHLH145 


Q9FGB0 


DNA binding 


Arabidopsis 
thaliana 


35.2/5.1 


18.3/5.1 


62 


28 


10/63 


92 


DNA-directed RNA 
polymerase subunit beta 


B1VKH5 


DNA-directed RNA 
polymerase activity 


Cryptomeria 
japonica 


139.6/9.4 


44.6/8.1 


59 


12 


11/65 


136 


DNA polymerase 
Protein synthesis 


P 10582 


3-5' exonuclease 
activity 


Zea mays 


108.2/8.6 


55.3/8.5 


58 


15 


8/23 


154 


50S ribosomal protein L33 

Protein folding, 
degradation and 
assembly 


B2LML4 


structural 
constituent of 
ribosome 


Guizotia 
abyssinica 


7.8/9.7 


56.4/5.1 


62 


42 


5/11 


45 


Heat shock protein 81-3 


P51818 


ATP binding 


Arabidopsis 
thaliana 


80.2/5.0 


26.7/8.2 


59 


19 


9/85 


73 


Proteasome subunit alpha 
type-6 

Transport 


048551 


threonine-type 
endopeptidase 
activity 


Glycine max 


27.5/5.8 


28.3/6.1 


58 


31 


6/62 


23 


NAD(P)H-quinone 

UXIUUlcUULldbc bUUUIIIl n, 

chloroplastic OS 


A6MMH1 


oxidoreductase 

dLUVIly, dLLIliy Ull 

NADH or NADPH 


Ch bra nth us 

C /Vv7 t\ I C 

bpiLULUb 


45.8/5.4 


18.5/7.3 


75 


26 


12/66 


JZ 


ndb leldLcU piULclil 

Mress Keiatea 


popi 07 
r Zo 1 0/ 


\j\r Ull lUll iy 


niauiuopsis 
thaliana 


Z^+. 1 / J.U 


Zo.D/o.^- 


O I 


ZO 


1 n/1 m 

I U/ I U I 


62 


Putative F-box/kelch- 
repeat protein At4g 19330 


Uco/U4 


M / A 


Arabidopsis 
thaliana 


62.9/7.1 


27.3/5.7 


59 


25 


1 0/82 


151 


Putative F-box/kelch- 
repeat protein At4g 19330 

Development 


065704 


N/A 


Arabidopsis 
thaliana 


62.8/7.1 


44.3/4.8 


65 


21 


8/75 


53 


3-ketoacyl-CoA synthase 5 


Q9C6L5 


fatty acid elongase 
activity 


Arabidopsis 
thaliana 


56.3/9.0 


27.7/8.1 


61 


40 


12/19 


81 


Cytoplasmic dynein 2 


Q9SMH5 


ATP binding 


Chlamydomonas 


483.5/6.1 


40.0/8.6 


58 


7 


22/48 




heavy chain 1 






reinhardtii 














Unknown 


















8 


Embryonic abundant 


P21745 


N/A 


Vicia faba 


30.1/6.4 


18.1/4.2 


67 


33 


7/67 




protein VF30.1 OS = Vicia 




















faba PE = 2 SV = 1 



















a As indicated in Figure 2. b Swissprot database accession number. c Molecular functions were inferred from those reported in the UniProt database. ^Theoretical 
isoelectric point and molecular mass. Experimental isoelectric point and molecular mass. ^Mascot score. ^Sequence coverage. ''Number of matched and 
unmatched peptides. 



could not be determined for 3% of the identified 1-1 pro- 
teins and -13% of the 1-2 proteins. The percentage of 
unidentified 1-2 proteins is similar to those found for 
maize [20], rice [21] or oilseed rape [22]. The substantial 
differences in the numbers of 1-1 and 1-2 proteins 
involved in metabolism, energy production, and protein 
synthesis suggest that the functions of these proteins dur- 
ing development deserve further attention. 

Subcellular localization provides important information 
about a protein's physiological function [23,24] . Recently, 



GO annotation has been widely used to predict the loca- 
tions of proteins [25-29], because the two are strongly cor- 
related. We found most of the identified proteins to be 
located in chloroplasts and mitochondria, which is congru- 
ent with total genomic data available for plants [30]. Inter- 
estingly, most of these proteins are involved in metabolism 
or energy production, which is consistent with the large 
number of proteins that we classified as metabolic or 
photosynthetic. We could not identify the cellular location 
of -11% and -25% 1-1 and 1-2 proteins, respectively. 
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Table 2 Identified proteins from the strawberry stolon 1-2 internode 



Spot Protein 
No. a 



Accession Molecular Function c 
Number 6 



Reference 
Organism 



Theoretical 
kDa/pl d 



Experimental 
kDa/pl e 



Score f SC 9 Matched/ 
(%) Unmatched 
queries' 7 



91 
95 

129 

139 

141 

172 

173 

174 



176 



17 



33 



35 



52 



64 



82 



90 



103 



104 



47 



121 



Metabolism 

Fructokinase-1 

Caffeic acid 3-0- 
methyltransferase 

Malate dehydrogenase, 
mitochondrial 

Glyceraldehyde-3-phosphate 
dehydrogenase, cytosolic 

Malate dehydrogenase, 
mitochondrial 

S-adenosylmethionine 
synthase 2 

S-adenosylmethionine 
synthase 

S-adenosylmethionine 
synthase 

Energy 

ATP synthase subunit beta, 
plastic 1 

Photosynthesis 

Probable alpha,alpha- 
trehalose-phosphate 
synthase [UDP-forming] 40S 

Ribulose bisphosphate 
carboxylase large chain 
(Fragment) OS 

Probable granule-bound 
starch synthase 1, 
chloroplastic/amyloplastic 

Ribulose bisphosphate 
carboxylase large chain 
(Fragment) 

Ribulose-phosphate 3- 
epimerase, chloroplastic 

Triosephosphate isomerase, 
chloroplastic 

Oxygen-evolving enhancer 
protein 1, chloroplastic 

Coproporphyrinogen-lll 
oxidase, chloroplastic 

Glutamyl-tRNA reductase 1, 
chloroplastic 

Transcription 

Homeobox-leucine zipper 
protein GLABRA 2 OS 

Protein HIRA 

Two-component response 
regulator ARR9 



128 DNA-directed RNA 

polymerase subunit beta 

149 B3 domain-containing 
protein Os07g0679700 

168 Protein HIRA 



A2WXV8 ATP binding 

Q8GU25 caffeate 0- 

methyltransferase activity 

P83373 Oxidoreductase 

P25858 Oxidoreductase 

P83373 Oxidoreductase 

Q9FUZ1 ATP binding 

A4PU48 ATP binding 

Q8W3Y4 ATP binding 



Q43157 

Q9M4S8 

P26320 

P35055 

Q42843 

P46607 
Q32SG6 
080366 

P 12465 

Q6Z3U3 

Q32SG6 



Oryza sotivo 
Rosa chinensis 



34.9/5.1 
40.1/5.6 



Frogorio 
ananasso 

Arabidopsis 
thaliana 

Frogario 
ononosso 



Medicogo 
truncatula 



35.8/8.7 
37.0/6.6 
35.8/8.7 
Brassico junceo 43.3/5.3 
43.7/5.6 
43.5/5.6 



Phaseolus 
lunatus 



Q8MBG5 ATP binding 



Cu scuta 
pentagona 



Q9T079 alpha,alpha-trehalose- Arabidopsis 

phosphate synthase thaliana 
(UDP-forming) activity 

P28261 magnesium ion binding Nypa fruticans 



Q9MAQ0 starch synthase activity 



P28391 magnesium ion binding 



ribulose-phosphate 3- 
epimerase activity 

triose-phosphate 
isomerase activity 

calcium ion binding 

coproporphyrinogen 
oxidase activity 

NADP or NADPH binding 



Arabidopsis 
thaliana 

Ceratopetalum 
gummiferum 

Spinacia 
oleracea 

Frag aria 
ananassa 

Solanum 
tuberosum 

Glycine max 



Hordeum 
vulgare 



sequence-specific DNA 
binding 

transcription regulator 
activity 

two-component 
response regulator 
activity 

DNA binding 

DNA binding 

transcription regulator 
activity 



Arabidopsis 
thaliana 

Zea mays 



Arabidopsis 
thaliana 



Chlorella 
vulgaris 

Oryza sativa 
Zea mays 



53.2/5.5 

90.3/6.1 

51.5/6.2 

67.5/8.7 

51.3/6.2 

30.6/8.2 
33.7/7.6 
35.6/5.8 
43.6/6.7 
58.1/8.7 

83.7/6.1 

106.55/7.8 

26.2/5.2 

1 79.4/9.9 
105.4/6.5 
106.5/7.8 



35.8/3.4 76 20 6/59 

43.7/4.4 110 35 12/30 

39.1/7.8 74 31 8/28 

40.1/9.3 74 29 9/27 

41.0/9.1 65 27 6/23 

43.9/5.9 138 44 15/36 

38.5/6.8 178 52 15/17 

41.2/6.7 132 48 15/40 

53.4/6.7 61 32 10/53 

14.2/6.3 58 25 12/66 

22.3/5.1 123 23 13/61 

23.0/4.5 60 24 9/37 

24.0/6.9 104 22 11/63 

27.2/7.9 58 31 6/51 

31.5/5.2 154 51 16/71 

34.5/4.7 98 31 10/48 

42.2/5.2 84 31 10/32 

42.1/4.9 70 24 11/48 

15.0/6.4 63 31 20/90 

19.7/7.1 60 22 16/90 

34.6/8.1 58 35 5/43 

36.9/8.0 60 9 11/34 

53.5/7.9 63 11 9/27 

41.1/6.8 60 22 16/90 
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Table 2 Identified proteins from the strawberry stolon 1-2 internode (Continued) 



Protein synthesis 

28 Eukaryotic translation 

initiation factor 5A-2 OS 

37 Eukaryotic translation 
initiation factor 5A 

45 Glutathione gamma- 

glutamylcysteinyltransferase 
2 

66 Phospho-2-dehydro-3- 

deoxyheptonate aldolase 1, 
chloroplastic 

87 Molybdenum cofactor 
sulfurase 



Q945F4 translation initiation 
factor activity 

Q9AXQ7 translation initiation 
factor activity 

Q9ZWB7 acyltransferase activity 



P21357 

Q655R6 
Q8RX55 



3-deoxy-7- 
phosphoheptulonate 
synthase activity 

lyase activity 



136 Probable beta-1, 3- Q8RX55 galactosyltransferase 

galactosyltransferase 18 activity 

Protein folding, 
degradation and assembly 

148 Anaphase-promoting Q8H1U5 
complex subunit 2 

153 U-box domain-containing Q8S8S7 
protein 34 

171 Ubiquitin carboxyl-terminal Q949Y0 
hydrolase 6 

175 U-box domain-containing Q8S8S7 

protein 34 activity 

Transport 

24 Magnesium transporter Q8H1G7 metal ion 
MRS2-8 OS 



Medicogo 
sotiva 

Di a nth us 
coryophyllus 

Arabidopsis 
thaliana 

Solon urn 
tuberosum 

Oryzo sativa 

Arabidopsis 
thaliana 



ubiquitin protein ligase Arabidopsis 

binding thaliana 

ubiquitin-protein ligase Arabidopsis 

activity thaliana 

ubiquitin-specific Arabidopsis 

protease activity thaliana 

ubiquitin-protein ligase Arabidopsis 
thaliana 



54 

4 
9 

20 
41 



Putative copper-transporting Q9SH30 
ATPase 3 

Stress Related 

Ninja-family protein AFP4 Q9S7Z2 
OS 

Glutathione S-transferase 6 Q96266 



Annexin D6 OS 



Monodehydroascorbate 
reductase, chloroplastic 

154 Monodehydroascorbate 
reductase 

Development 

2 Pentatricopeptide repeat- 
containing protein 
At3g06430, chloroplastic OS 

Protein BRUSHY 1 



Q9LX08 
P92947 
Q40977 



transmembrane 
transporter activity 

ATP binding 



protein binding 
glutathione binding 
calcium ion binding 



Arabidopsis 
thaliana 



Arabidopsis 
thaliana 



Arabidopsis 
thaliana 

Arabidopsis 
thaliana 

Arabidopsis 
thaliana 



17.5/5.4 
17.6/5.6 
52.3/6.6 

60.0/8.9 

92.9/7.1 
77.9/8.7 

98.4/4.8 
91.7/9.1 
54.0/5.8 
91.7/9.1 

43.1/5.3 
109.0/6.0 

35.6/8.5 
29.2/8.5 
36.6/7.7 
53.5/8.1 



monodehydroascorbate Arabidopsis 

reductase (NADH) activity thaliana 

monodehydroascorbate Pisum sativum 47.3/5.8 
reductase (NADH) activity 



11 

38 



12 
16 



Q9SQU6 N/A 



Q6Q4D0 protein binding 



Arabidopsis 
thaliana 



Arabidopsis 
thaliana 



Protein PAIR1 

1 -aminocyclopropane-1 - 
carboxylate synthase 7 



Q75RY2 
Q9STR4 



Oryza sativa 

1 -aminocyclopropane-1 - Arabidopsis 
carboxylate synthase thaliana 
activity 



Unknown 

Probable protein ABIL5 OS Q5JMF2 

Pentatricopeptide repeat- Q9SZ20 
containing protein 
At4g26800 OS 

Putative F-box protein Q9LM74 
At1g20795 OS 



N/A 
N/A 

N/A 



Oryza sativa 

Arabidopsis 
thaliana 

Arabidopsis 
thaliana 



56.2/7.8 

148.8/5.5 

53.8/9.8 
51.0/5.9 



28.3/8.3 
58.4/9.1 

48.5/8.6 



1 7.3/4.7 
20.3/5.9 
16.4/6.6 

27.9/7.3 

27.2/5.2 
35.1/9.1 

43.5/7.9 
39.6/6.9 
41.6/5.8 
44.7/6.2 

13.6/4.4 
29.0/8.0 

7.1/5.0 
1 1 .9/6.2 
14.0/5.1 
16.2/5.7 
43.5/6.8 

9.5/6.7 

12.0/5.9 

1 7.2/5.5 
26.2/6.3 

15.2/7.6 
16.1/6.8 

20.3/5.1 



74 
88 
57 

62 

60 
66 

64 
63 
58 
58 



57 
58 
62 
63 
58 



63 

66 
57 



26 5/39 

45 7/28 

15 8/32 

31 11/71 

12 6/22 

1 7 7/27 

12 7/20 

18 12/51 

19 8/35 
19 12/67 



58 22 6/34 



65 18 14/48 



40 9/70 

25 5/46 
20 7/39 

26 8/76 
22 7/31 



58 25 11/59 



24 21/71 



19 10/34 
19 10/72 



59 40 7/42 
64 32 12/49 

60 21 8/55 
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Table 2 Identified proteins from the strawberry stolon 1-2 internode (Continued) 



67 



113 



Probable protein ABIL5 OS Q5JMF2 N/A 

BRCT domain-containing 004251 N/A 
protein At4g021 10 

Thylakoid lumenal 15.0 kDa Q9LW5 N/A 
protein 2, chloroplastic 

Pentatricopeptide repeat- Q5XET4 N/A 

containing protein 

At2g01860 



Oryza sotivo 

Arobidopsis 
thaliana 

Arabidopsis 
thaliana 

Arabidopsis 
thaliana 



28.2/8.2 
142.5/8.4 

24.7/5.7 

56.0/9.2 



21.5/6.9 
28.4/7.5 

37.8/5.8 

36.1/6.8 



59 
59 



59 



40 7/42 

19 18/80 

26 7/37 

39 15/95 



a Spot number as in Figure 2. b Swissprot database accession number. c Molecular functions were inferred from those reported in the UniProt database. ^Theoretical 
molecular mass and isoelectric point. Experimental molecular mass and isoelectric point. ^Mascot score. ^Sequence coverage. ''Number of matched and unmatched 
peptides. 



Proteins involved in metabolism 

For 1-1, identified proteins were found for lipid metabo- 
lism (spot 42), purine metabolism (spot 40), the citric 
acid cycle (spots 15, 90, and 91), starch and sucrose 
metabolism (spots 29, 31, 118, 133, and 162), and one- 
carbon metabolism (spots 111 and 142). For 1-2, we 
found additional proteins involved in one-carbon meta- 
bolism (spots 172, 173, and 174) but fewer and different 
proteins involved in starch and sucrose metabolism 
(spots 91 and 139). Many more proteins involved in 
metabolism were found for 1-1 than for 1-2 (Tables 1 
and 2, Figure 5). The activities of extracellular lipase 
(spot 42) and P-l,3-glucanase (spots 29 and 31) have 
been implicated in pollen germination, fertilization, 
response to wounding, and cell division [31,32], all or 
any of which may be related to 1-1 elongation via cell 
division. Sugars are involved in energy metabolism and 
act as signaling molecules. For 1-1, we identified 2,3- 
bisphosphoglycerate-independent phosphoglycerate 
mutase (spot 118) and phosphoglucomutase (spot 162), 
and for 1-2, glyceraldehyde-3-phosphate dehydrogenase 
(spot 159), all of which are glycolytic enzymes. Interest- 
ingly, malate dehydrogenase (spots 90 and 129), fructoki- 
nase-1 (spots 155 and 91), and S-adenosyl-l-methionine 




Figure 3 Venn diagram for the identified 1-1 and 1-2 proteins. 

The numbers and percentages of unique proteins (excluding 
isoforms) found for either or both internodes are given. Spot 
numbers for 1-1 proteins are given first followed by spot numbers 
for the corresponding I-2 proteins in parentheses. 



(AdoMet) synthase (spots 111 and 175) were identified in 
both 1-1 and 1-2 proteomes. Malate dehydrogenase is a 
citric-acid-cycle enzyme that converts malate into oxa- 
loacetate (using NAD). Fructokinase is involved in 
sucrose and fructose metabolism and may regulate starch 
synthesis in conjunction with sucrose synthase, which 
first metabolizes plant sink tissue in, for example, pota- 
toes [33]. AdoMet synthase catalyzes the formation of 
AdoMet from methionine and ATP [34], which is the 
main methyl group donor and is involved in transmethy- 
lations and the trans-sulfuration pathway [35]. AdoMet is 
also involved in the biosynthesis of many secondary 
metabolites [36,37] and can be decarboxylated to gener- 
ate a propylamine donor for polyamine biosynthesis [38]. 
Polyamines are required for cell proliferation and may 
play a role in the rapid growth of bloom-forming dinofla- 
gellates [39]. In plants, AdoMet participates in ethylene 
biosynthesis [40] and is the methyl group donor in trans- 
methylation of alkaloids [41]. Cell and life cycle variation 
in AdoMet synthase expression has been observed in 
yeast and apicomplexa [42,43]. 

Proteins involved in energy production 

More identified 1-1 proteins (11%) were found to be 
involved in energy production than were 1-2 proteins 
(2%, Figure 5). Both the a and P subunits of chloroplastic 




25 

I 20 

4) 

£ 15 

8 10 
1 

* 5 



Arabidopsis Oryza sativa 
thaliana 

Figure 4 Number of proteins with sequences that matched 
those of organisms listed in the Swissprot Viridiplantae 
database. Over 50% of the proteins identified had sequences 
similar to annotated proteins from dicot species in the Viridiplantae 
database. 
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□ Metabolism iEnergy □ Photosynthesis □ Transcription 

■ Protein synthesis DProtein folding ■ Transport □ Stress 

■ Development ■ Unknown 

Figure 5 Pie charts classifying the identified 1-1 and 1-2 proteins according to biological function The identified proteins were grouped 
according to their biological processes and are expressed in percentage. 



ATP synthase (spots 156 and 152) were identified for 1-1, 
whereas only the ATP synthase P-subunit was detected 
in the 1-2 proteome. ATP synthase is a very large com- 
plex (>500 kDa) embedded in the inner membranes of 
chloroplasts and mitochondria. It utilizes the products of 
fat and carbohydrate breakdown to generate proton gra- 
dients across membranes, which then drive ATP 
synthesis. 

Proteins involved in photosynthesis 

Photosynthesis uses light energy and chlorophyll to 
synthesize simple sugars from carbon dioxide and water 
and to capture the energy as phosphate bonds in ATP. 
ATP is then available as an energy source, and the 
sugars are used as building blocks to produce other cell 
structural and storage components. Photosynthetic 
enzymes including RuBisCo (spots 2, 37, 46, and 33), 
oxygen- evolving enhancer protein 1 (spots 126 and 90), 




Figure 6 Gene ontology classification of identified 1-1 and 1-2 
proteins according to their subcellular location. Subcellular 
locations of the proteins were assigned according to the GO 
annotations and are expressed as percentages of the assigned 
proteins. 



and triosephosphate isomerase (spots 69, 72, and 82) 
were found in the 1-1 and 1-2 proteomes. 

Interestingly, many other proteins involved in sugar 
synthesis were found only in the 1-2 proteome, i.e., 
alpha, alpha-trehalose-phosphate synthase (UDP-form- 
ing, spot 17), granule-bound starch synthase 1 (spot 
35), and ribulose-phosphate 3-epimerase (spot 64). 
Alpha, alpha-trehalose-phosphate synthase (UDP-form- 
ing) synthesizes alpha, alpha-trehalose 6-phosphate 
from D-glucose 6-phosphate, and is then dephosphory- 
lated to trehalose by trehalose 6-phosphate phospha- 
tase. Trehalose metabolism, a side-branch of carbon 
flux in bacteria, yeast, and plants, has recently drawn 
attention because it may partially regulate plant 
growth, development, and stress resistance [44]. Gran- 
ule-bound starch synthase 1 is required for the synth- 
esis of amylase, and ribulose-phosphate 3-epimerase 
(pentose-5-phosphate 3-epimerase) converts d-ribulose 
5-phosphate into d-xylulose 5-phosphate as part of the 
reductive pentose phosphate (Calvin) cycle. These 
aforementioned enzymes are very important for sugar 
synthesis. Sugars can act as signaling molecules in 
microorganisms, animals, and plants. During plant 
growth and development, sugars modulate seed germi- 
nation, seedling development, root and leaf differentia- 
tion, floral transition, fruit ripening, embryogenesis, 
senescence, and responses to light, stress, and patho- 
gens [45-53]. For strawberries, the clonal plant is 
usually found at the second node, so the identification 
of specific sugar-related enzymes found only in the 1-2 
proteome should increase our understanding of the 
different internodal developmental mechanisms. Our 
results indicate that positive interactions between 
sugar synthesis and hormonal signaling in 1-2 may be 
necessary for asexual strawberry reproduction. 
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Proteins involved in protein synthesis 

We identified only one 1-1 protein associated with pro- 
tein synthesis, namely the 50S ribosomal protein L33 
(spot 154). Cell growth and division require the synth- 
esis of new proteins and ribosomes. L33 is involved in 
the biogenesis of both the small and large ribosomal 
subunits. For 1-2, we found the eukaryotic translation 
initiation factor 5 A (spots 28 and 37), which is a highly 
conserved eukaryotic protein. Eukaryotic translation 
initiation factor 5A appears to be involved in RNA 
metabolism and trafficking in mammals and yeast, thus 
regulating cell proliferation, cell growth, and pro- 
grammed cell death [54]. In plants, however, its physio- 
logical function is not known. 

Proteins involved in protein folding and protein 
degradation 

For 1-1, two proteins involved in protein folding and 
processing were identified, and for 1-2, four proteins 
involved in protein modification and degradation were 
found. The heat shock protein 81-3 (spot 45) found in 
the 1-1 proteome is a molecular chaperone and likely 
involved in signal transduction and development asso- 
ciated with certain hormone receptors and kinases [55]. 
Interestingly, three of the 1-2 proteins, anaphase-pro- 
moting complex subunit 2 (APC2; spot 148), U-box 
domain-containing protein 34 (spots 153 and 175), and 
ubiquitin carboxyl-terminal hydrolase 6 (spot 171), are 
all ubiquitin-dependent proteins involved in catabolism. 
The ubiquitin-proteasome pathway is responsible in 
large part for protein degradation and consequently reg- 
ulates many aspects of development. The identification 
of the aforementioned proteins suggests that research 
on the ubiquitin conjugation pathway might illuminate 
the mechanism of 1-2 clonal multiplication. 

Proteins involved in transcription 

We identified five 1-1 proteins, histone acetyltransferase 
GCN5 (spot 6), transcription factor bHLH145 (spot 27), 
DNA polymerase (spot 136), pentatricopeptide repeat- 
containing protein (spot 18), and DNA-directed RNA 
polymerase (spot 92), that are involved in transcription. 
GCN5 is a coactivator of transcriptional regulation 
[56,57]. The 1-2 homeodomain-leucine zipper protein 
(spot 18) is a putative transcription factor required for 
correct morphological development and maturation of 
trichomes as well as for normal development of seed 
coat mucilage [58]. The function of histone regulator 
protein (spots 47 and 168) has yet to be determined; 
however, it may be involved in maintining knox genes 
silencing throughout leaf development [59]. The two- 
component response regulator (spot 121) is involved in 
the His-to-Asp phosphorelay signal transduction system 
[60], 



Proteins involved in stress response and development 

Twelve of the identified 1-1 and 1-2 proteins are asso- 
ciated with stress and development. One 1-1 protein 
(spots 62 and 151) and five 1-2 proteins (e.g., spots 4, 9, 
20, 41, and 154) are involved in stress responses, and 
two 1-1 proteins (spots 53 and 81) and four 1-2 proteins 
(spots 2, 11, 38, and 68) are associated with develop- 
ment. Spot 53 is 3-ketoacyl-CoA synthase, which med- 
iates the synthesis of very long chain fatty acids (26 to 
30 carbons). Spot 81 is the cytoplasmic dynein 2 heavy 
chainl, which is an intracellular motor for retrograde 
vesicle and organelle motility along microtubules. 

For 1-2, monodehydroascorbate reductase (spots 41 
and 154) was identified and is an oxidoreductase that 
oxidizes NADH or NADPH using a quinone as the oxi- 
dant during the glutathione-ascorbate cycle, a major 
plant antioxidant system that protects against reactive 
oxygen species. Monodehydroascorbate reductase activ- 
ity has been found in chloroplasts, the cytosol, mito- 
chondria, glyoxysomes, and leaf peroxisomes [61]. 1- 
aminocyclopropane-l-carboxylate synthase 7 (spot 68) 
catalyzes the conversion of AdoMet into 1-aminocyclo- 
propane-l-carboxylate, a precursor of ethylene. 

Comparison of metabolic pathways in 1-1 and 1-2 

Interestingly, many of the proteins involved in central 
metabolic pathways (e.g., glycolysis, the citric acid cycle, 
pyruvate metabolism) were identified in both the 1-1 
and 1-2 proteomes. Using the KEGG PATHWAY data- 
base (http://www.genome.jp/kegg/pathway.html), we 
classified more of the 1-1 proteins than the 1-2 proteins 
as involved in carbon fixation, glyoxylate and dicarboxy- 
late metabolism, glycolysis/gluconeogenesis, and oxida- 
tive phosphorylation (Figure 7). For the two internodes, 
seven enzymes (15 spots) were classified as carbon fix- 
ing, and 13.2% of the 1-1 proteins belonged to this cate- 
gory, whereas only 7.7% of those from 1-2 did. The 
difference is related to the number of spots found for 
the RuBisCo complex, i.e., more were found for 1-1 than 
for 1-2. Most of the enzymes of the citrate cycle, pyru- 
vate metabolism, fructose and mannose metabolism, and 
starch and sucrose metabolism were identified in the 
two proteomes, and the numbers of proteins found for 
each pathway were similar, suggesting that the house- 
keeping pathways are needed for stolon viability. Many 
enzymes involved in the ubiquitin-proteasome pathway 
were also identified but were greater in relative number 
in the 1-2 proteome (Figure 7). In plants, regulated pro- 
tein degradation by the ubiquitin-proteasome system 
contributes substantially to development by affecting 
many processes, e.g., embryogenesis, hormone signaling, 
and senescence, which suggests that the ubiquitin-pro- 
teasome system may play a central role in 1-2 morpho- 
genesis. The 1-1 and 1-2 proteomes, which we have 
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Protein abundance (t2-DE gel voline) 

Figure 7 Protein abundance for nine metabolic pathways in 1-1 

and 1-2. For carbon fixation, glyoxylate and dicarboxylate 
metabolism, glycolysis/gluconeogenesis, pentose phosphate 
pathway, pyruvate metabolism, starch and sucrose metabolism, 
oxidative phosphorylation, ubiquitin-proteasome pathway; 7, 6, 4, 1, 
1, 3, 6, and 5 enzymes were identified, respectively, corresponding 
to 15, 10, 6, 5, 5, 4, 6, and 6 2-DE spots. 

V J 



described herein, will allow us to analyze and compare 
proteins in metabolic pathways and may provide new 
insights into the regulation and expression of various 
molecular networks. 

Possible PPI networks and adventitious root-formation 
mechanisms in 1-2 

Proteins are the main catalysts, structural elements, sig- 
naling messengers, and molecular machines of biological 




Figure 8 Possible protein-protein interaction network among I- 
2 proteins derived using the Cytoprophet module of 
Cytoscape. Cytoprophet draws a network of potential interactions 
with probability scores and GO distances as edge attributes. 
Proteins are marked with UniProt ID names. 

v / 



tissues [62]. PPIs are extremely important in orchestrat- 
ing cellular events. Protein interaction networks provide 
road maps of cellular pathways. Therefore, many meth- 
ods are used to characterize PPIs, which include physi- 
cal interactions and functional linkages [63]. 

Although we could not use Cytoprophet to delineate a 
PPI network for 1-1 because of limited information, a pos- 
sible PPI network was found for 1-2 (Figure 8). U-box 
domain-containing protein 34 (PUB34) (Q8S8S7_ARATH) 
is the central core protein of the signaling network, as it 
interacts with many other proteins, e.g., alpha, alpha-tre- 
halose-phosphate synthase (TPS4) (Q9T079_ARATH), 
APC2 ( Q8 H 1 U5_ AR ATH) , 1-aminocyclopropane-l-car- 
boxylate synthase 7 (ACS7) (1A17_ ARATH), and EIF5A. 
APC2 is an E3 ubiquitin ligase that is a component of the 
SCF family ubiquitin ligases, which catalyze the attach- 
ment of ubiquitin to the lysine side chains of securin and 
mitotic cyclins [64-66]. EIF5A interacts with BRUSHY1 
(BRU1_ARATH), which is required for the proper 
arrangement of cells in the root and shoot apical meris- 
tems. Ubiquitin-mediated protein degradation probably 
affects meristem structural formation by modulating the 
concentration of cell-cycle regulators and transcription 
factors [67]. Therefore, the ubiquitin system may be vital 
during the morphogenesis of clonal plants. 

Adventitious roots develop from the second node at 
the end of 1-2 before clonal plant formation. We devel- 
oped a model for adventitious root formation in 1-2 
based on published data [64-69] and our findings; the 
model includes four regulated pathways (Figure 9). 
Regulated protein degradation has repeatedly been iden- 
tified as a key component of cell-cycle regulation. 
Securin inhibits a protease called separase, which cleaves 
cohesins allowing anaphase onset. Activated APC cdc20 
targets securin for degradation, which initiates the meta- 
phase-to-anaphase transition [68]. In addition, biochem- 
ical and molecular studies have shown that EIF-5A is 
crucial for plant growth and development as it regulates 
cell division and cell growth [69]. Continuous cell divi- 
sion, elongation, and differentiation can cause the for- 
mation of root primordia, so the APC complex-related 
and EIF-5A-related biological processes may be two 
important pathways that regulate the formation of 
adventitious roots. Moreover, ACS7 catalyzes the con- 
version of AdoMet into 1-aminocyclopropane-l-carbox- 
ylate, a direct precursor of ethylene, whereas ACS7 is 
ubiquitinated. Ubiquitination probably leads to its subse- 
quent degradation, thus controlling ethylene production. 
Ethylene can regulate root initiation and emergence. 
Conversely, as an important catalytic enzyme, TPS4 
plays a central role in the complex signaling network 
that links sugars and hormones with its interacting part- 
ner PUB34. Together, the four pathways work synergis- 
tically to induce formation of adventitious roots. 
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Adventitious root formation of a new clonal plant 



Figure 9 Model for adventitious root and clonal plant 
formation in 1-2 that incorporates four regulated pathways. 

Five identified 1-2 proteins were integrated into the model, and the 

possible PPIs are shown (dashed lines) based on the PPI network in 

Figure 8. (a) Anaphase-promoting complex (APC/C) is a ubiquitin 

ligase that plays a key role in the cell cycle, (b) Eukaryotic translation 

initiation factor 5A (EIF-5A) may interact with PUB34 to regulate cell 

division, (c) ACS7, when interacting with ubiquitin ligase, plays a 

central role in ethylene biosynthesis, (d) Important regulatory effects 

on plant growth and development have been reported for 

trehalose (Tre) and trehalose 6-phosphate (T6P). (CDC20: cell-division 

cycle protein 20; G6P: glucose 6-phosphate; TPP: trehalose-6- 

phosphate phosphatase 
v ) 



Conclusions 

For this research, we compared the proteomes of the 1-1 
and 1-2 internodes of the strawberry stolon to begin to 
elucidate the how the differences in the proteomes affect 
the growth and functional characteristics of the two 
internodes as only the second node, at the end of 1-2, 
tends to form adventitious roots and a clonal plant. In I- 
1, the majority of the proteins were involved in metabo- 
lism, photosynthesis, energy, and transcription. In 1-2, 
relatively more proteins were involved in photosynthesis, 
carbohydrate metabolism, stress responses, and protein 
folding and degradation, indicating that these many dif- 
ferent processes work synergistically to induce cell differ- 
entiation necessary for root and plant formation. Given 
our findings and those of others, we present a scheme for 
protein interactions that could be responsible for adven- 
titious root and clonal plant formation in 1-2. 

Acknowledgements 

This research was supported in part by the Key Scientific and Technological 
Innovation Projects of Hangzhou (Grant No. 2007231 2A03 to Songlin Ruan). 

Author details 

laboratory of Plant Molecular Biology and Proteomics, Institute of Biology, 
Hangzhou Academy of Agricultural Sciences, Hangzhou 310024, China. 



2 School of Life Science, Zhejiang Chinese Medical University, Hangzhou 
310053, China. 

Authors' contributions 

XF designed and carried out experiments, analyzed data, and wrote the 
manuscript. HM designed the study, interpreted data, and assisted with the 
writing of the manuscript. DL did most of the MS analysis and assisted with 
the evaluation of MS data. HY conceived the study and participated in its 
design. WL designed and analyzed experiments. SR conceived of, designed, 
and coordinated the study, and assisted with the writing of the manuscript. 
All authors read and approved the final manuscript. 

Competing interests 

The authors declare that they have no competing interests. 

Received: 21 January 2011 Accepted: 14 May 2011 
Published: 14 May 2011 

References 

1. Stuefer JF: Separating the effects of assimilate and water integration in 
clonal fragments by the use of steam-girdling. Abstr Bot 1995, 19:75-81. 

2. Stuefer JF, Huber H: The role of stolon internodes for ramet survival after 
clone fragmentation in potentilla anserina. Ecol Lett 1999, 2:135-139. 

3. Alonso-Blanco C, Mendez-Vigo B, Koornneef M: From phenotypic to 
molecular polymorphisms involved in naturally occurring variation of 
plant development. Int J Dev Biol 2005, 49:71 7-732. 

4. Jomn JD, Rubiales E, Dumas-Gaudot AM, Maldonado M, Castillejo A, 
Curto M: Proteomics: a promising approach to study biotic interaction in 
legumes. A review. Euphytico 2006, 147:37-47. 

5. Bianco L, Lopez L, Scalone AG, Carli MD, Desiderio A, Benvenuto E, 
Perrotta G: Strawberry proteome characterization and its regulation 
during fruit ripening and in different genotypes. J Proteomics 2009, 
72:586-607. 

6. Giavalisco P, Nordhoff E, Kreitler T, Kloppel KD, Lehrach H, Klose J, Gobom J: 
Proteome analysis of Arabidopsis thaliana by two-dimensional gel 
electrophoresis and matrix-assisted laser desorption/ionisation-time of 
flight mass spectrometry. Proteomics 2005, 5:1902-1913. 

7. Larrainzar E, Wienkoop S, Wolfram W, Ruben L, Cesar Al, Gonzalez EM: 
Medicago truncatula root nodule proteome analysis reveals differential 
plant and bacteroid responses to drought stress. Plant Physiol 2007, 
144:1495-1507. 

8. Hjerno K, Aim R, Canback B, Matthiesen R, Trajkovski K, Bjork L, Roepstorff P, 
Emanuelsson CS: Down-regulation of the strawberry Bet v 1 -homologous 
allergen in concert with the flavonoid biosynthesis pathway in 
colourless strawberry mutant. Proteomics 2006, 6(Suppl 5):1 574-1 587. 

9. Takashi N, Yutaka H: Elongation of strawberry runners in relation to 
length and number of cells. J Agricultural Res 1993, 43:3-4. 

10. Bradford MM: A rapid and sensitive method for quantitation of 
microgram quantities of protein utilizing the principle of protein-dye 
binding. Anal Biochem 1976, 72:248-254. 

11. Wheeler CH, Dunn MJ, Rheumatology H, London UK, Pharmacia A: A 
modified silver staining protocol for visualization of proteins compatible 
with matrix-assisted laser desorption/ionization and electrospray 
ionization-mass spectrometry. Electrophoresis 2000, 21:3666-3672. 

12. Mechin V, Balliau T, Chateau-Joubert S, Davanture M: A two-dimensional 
proteome map of maize endosperm. Phytochemistry 2004, 65:1609-1618. 

13. Giavalisco P, Nordhoff E, Kreitler T, Kloppel KD: Proteome analysis of 
Arabidopsis thaliana by two-dimensional gel electrophoresis and matrix- 
assisted laser desorption/ionisation-time of flight mass spectrometry. 
Proteomics 2005, 5:1902-1913. 

14. Jiang Y, Yang B, Harris NS, Deyholos MK: Comparative proteomic analysis 
of NaCI stress-responsive proteins in Arabidopsis roots. J Exp Bot 2007, 
58:3591-3607. 

15. Ishikawa T, Yoshimura K, Tamoi M, Takeda T, Shigeoka S: Alternative mRNA 
splicing of 3'-terminal exons generates ascorbate peroxidase isoenzymes 
in spinach (Spinacia oleracea) chloroplasts. J Biochem 1997, 328(Suppl 
3):795-800. 

16. Holmes-Davis R, Tanaka CK, Vensel WH, Hurkman WJ, McCormick S: 
Proteome mapping of mature pollen of Arabidopsis thaliana. Proteomics 
2005, 5(Suppl 18)4864-4884. 



Fang et al. Proteome Science 201 1, 9:26 
http://www.proteomesci.eom/content/9/1/26 



Page 14 of 15 



17. Ndimba BK, Chivasa S, Simon WJ, Slabas AR: Identification of Arabidopsis 
salt and osmotic stress responsive proteins using two-dimensional 
difference gel electrophoresis and mass spectrometry. Proteomics 2005, 
5(Suppl 16):41 85-41 96. 

18. Huber SC, Hardin SC: Numerous post translational modifications provide 
opportunities for the intricate regulation of metabolic enzymes at 
multiple levels. Curr opin Plant Biol 2004, 7:318-322. 

19. Hoffmann R, Valencia A: A gene network for navigating the literature. Nat 
Genet 2004, 36:664. 

20. Mechin V, Balliau T, Chateau-Joubert S, Davanture M, Langella 0, Negroni L, 
Prioul JL, Thevenot C, Zivy M, Damerval M: A two-dimensional proteome 
map of maize endosperm. Phytochemistry 2004, 65:1609-1618. 

21. Koller A, Washburn MP, Lange BM, Andon NL, Deciu C, Haynes PA, Hays L, 
Schieltz D, Ulaszek R, Wei J, Wolters D, John R, Yates JR: Proteomic survey 
of metabolic pathways in rice. Proc Natl Acad Sci USA 2002, 
18:11969-11974. 

22. Hajduch M, Casteel JE, Hurrelmeyer KE, Song Z, Agrawal GK, Thelen JJ: 
Proteomic analysis of seed filling in Brassica napus. developmental 
characterization of metabolic isozymes using high resolution two- 
dimensional gel electrophoresis. Plant Physiol 2006, 141:32-46. 

23. Cocco L, Manzoli L, Barnabei 0, Martelli AM: Significance of subnuclear 
localization of key players of inositol lipid cycle. Adv Enzyme Regul 2004, 
44:51-60. 

24. Itoh K, Brott BK, Bae GU, Ratcliffe MJ, Sokol SY: Nuclear localization is 
required for dishevelled function in wnt/beta-catenin signaling. J Biol 
2005, 4(Suppl 1):3. 

25. Chou K, Cai Y: Prediction of protein subcellular locations by GO-FunD- 
PseAA predictor. Biochem Biophys Res Commun 2004, 320:1236-1239. 

26. Lei Z, Dai Y: Assessing protein similarity with gene ontology and its use 
in subnuclear localization prediction. BMC Bioinf 2006, 7:491. 

27. Blum T, Briesemeister S, Kohlbacher O: MultiLoc2: integrating phylogeny 
and gene ontology terms improves subcellular protein localization 
prediction. BMC Bioinf 2009, 10:274. 

28. Briesemeister S, Blum T, Brady S: SherLoc2: a high-accuracy hybrid 
method for predicting subcellular localization of proteins. J Proteome Res 
2009, 8:5363-5366. 

29. Huang WL, Tung CW, Ho SW: ProLoc-GO: utilizing informative gene 
ontology terms for sequence-based prediction of protein subcellular 
localization. BMC Bioinf 2008, 9:80. 

30. Van Wijk KJ: Challenges and prospects of plant proteomics. Plant Physiol 
2001, 126:501-508. 

31. Brick DJ, Brumlik MJ, Buckley JT, Cao JX, Davies PC, Misra S, Tranbarger TJ, 
Upton C: A new family of lipolytic plant enzymes with members in rice, 
arabidopsis and maize. FEBS Lett 1995, 377:475-480. 

32. Abeles FB, Bosshart RP, Forrence LE, Habig WH: Preparation and 
purification of glucanase and chitinase from bean leaves. Plant Physiol 
1970, 47:129-134. 

33. Odanaka S, Bennett AB, Kanayama Y: Distinct physiological roles of 
fructokinase isozymes revealed by gene-specific suppression of Frk1 and 
Frk2 expression in tomato. Plant Physiol 2002, 20:703-706. 

34. Catoni GL: S-Adenosylmethionine: a new intermediate formed 
enzymatically from L-methionine and adenosinetriphosphate. J Biol 
Chem 1953, 204(suppl 1):403-416. 

35. Chiang PK, Gordon RK, Tal J, Zeng GC, Doctor BP: Pardhasaradhi K. 
McCann PP: S-Adenosylmethionine and methylation. FASEB J 1996, 
10(suppl 4):471-480. 

36. Shimizu Y, Gupta S, Masuda K, Maranda L, Walker CK, Wang R: 
Dinoflagellate and other microalgal toxins: chemistry and biochemistry. 
Pure Appl Chem 1989, 61(Suppl 3):513-516. 

37. Shimizu Y: Microalgal metabolites: a new perspective. Annu Rev Microbiol 
1996, 50:431-465. 

38. Tabor H, Rosenthal SM, Tabor CW: The biosynthesis of spermidine and 
spermine from putrescine and methionine. J Biol Chem 1958, 233(Suppl 
4):907-914. 

39. Chan KL, New D, Ghandhi S, Wong F, Lam CM, Wong JT: Transcript levels 
of the eukaryotic translation initiation factor 5A gene peak at early G(1) 
phase of the cell cycle in the dinoflagellate Crypthecodinium cohnii. 

Appl Environ Microbiol 2002, 68(Suppl 5):2278-2284. 

40. Yang S, Hoffman N: Ethylene biosynthesis and its regulation in higher 

plants. Annu Rev Plant Physiol 1984, 35:155-189. 



41. Kutchan TM: Alkaloid biosynthesis-the basis for metabolic engineering of 
medicinal plants. Plant Cell 1995, 7(Suppl 7):1 059-1 070. 

42. Slapeta J, Stejskal F, Keithly JS: Characterization of S-adenosylmethionine 
synthetase in Cryptosporidium parvum (Apicomplexa). FEMS Microbiol 
Lett 2003, 225(Suppl 2):271-277. 

43. Brejning J, Jespersen L: Protein expression during lag phase and growth 
initiation in Saccharomyces cerevisiae. Int J Food Microbiol 2002, 75(Suppl 
1-2):27-38. 

44. Rolland F, Winderickx J, Johan M: Glucose-sensing mechanisms in 
eukaryotic cells. Trends Biochem Sci 2001, 26:310-317. 

45. Koch KE: Carbohydrate modulated gene expression in plants. Annu Rev 
Plant Physiol Plant Mol Biol 1996, 47:509-540. 

46. Roitsch T: Source-sink regulation by sugar and stress. Curr Opin Plant Biol 
1999, 2:198-206. 

47. Sheen J, Zhou L, Jang JC: Sugars as signaling molecules. Curr Opin Plant 
Biol 1999, 2:410-418. 

48. Wobus U, Weber H: Sugars as signal molecules in plant seed 
development. Biol Chem 1999, 380:937-944. 

49. Smeekens S: Sugar-induced signal transduction in plants. Annu Rev Plant 
Physiol Plant Mol Biol 2000, 51:49-81. 

50. Xiao W, Sheen J, Jang JC: The role of hexokinase in plant sugar signal 
transduction and growth and development. Plant Mol Biol 2000, 
44:451-461. 

51. Coruzzi GM, Zhou L: Carbon and nitrogen sensing and signaling in 
plants: emerging matrix effects. Curr Opin Plant Biol 2001, 4:247-253. 

52. Ohto M, Onai K, Furukawa Y, Aoki E: Effects of sugar on vegetative 
development and floral transition in Arabidopsis. Plant Physiol 2001, 
127:252-261. 

53. Corbesier L, Bernier G, Perilleux C: C:N ratio increases in the phloem sap 
during floral transition of the long-day plants Sinapis alba and 
Arabidopsis thaliana. Plant Cell Physiol 2002, 43:684-88. 

54. Kang HA, Schwelberger HG, Hershey JW: Translation initiation factor elF- 
5A, the hypusine-containing protein, is phosphorylated on serine in 
Saccharomyces cerevisiae. J Biol Chem 1993, 268:14750-14756. 

55. Saracco SA, Hansson M, Scalf M, Walker JM, Smith LM, Vierstra RD: Tandem 
affinity purification and mass spectrometric analysis of ubiquitylated 
proteins in Arabidopsis. J Plant 2009, 59:344-358. 

56. Brownell JE, Zhou J, Ranalli T, Kobayashi R, Edmondson DG, Roth SY, 
Allis CD: Tetrahymena histone acetyltransferase A: a homolog to yeast 
Gcn5p linking histone acetylation to gene activation. Cell 1996, 84(Suppl 
6):843-851. 

57. Lee Tl, Young RA: Transcription of eukaryotic protein-coding genes. Annu 
Rev Genet 2000, 34:77-137. 

58. Ohashi Y, Oka A, Ruberti I, Morelli G, Aoyama T: Entopically additive 
expression of GLABRA2 alters the frequency and spacing of trichome 
initiation. J Plant 2002, 29:359-369. 

59. Phelps-Durr TL, Thomas J, Vahab P, Timmermans MCP: Maize rough 
sheath2 and its Arabidopsis orthologue ASYMMETRIC LEAVES 1 interact 
with HIRA, a predicted histone chaperone, to maintain knox gene 
silencing and determinacy during organogenesis. Plant Cell 2005, 
17:2886-2898. 

60. To JPC, Haberer G, Ferreira FJ, Deruere J, Mason MG, Schaller GE, Alonso JM, 
Ecker JR, Kieber JJ: Type-A Arabidopsis response regulators are partially 
redundant negative regulators of cytokinin signaling. Plant Cell 2004, 
16:658-671. 

61. Leterrier M, Corpas F, Barroso J, Sandalio L: Peroxisomal 
monodehydroascorbate reductase. Genomic clone characterization and 
functional analysis under environmental stress conditions. Plant Physiol 
2005, 138(Suppl 4):21 11-2123. 

62. Eisenberg D, Marcotte EM, Xenarios I, Yeates TO: Protein function in the 
post-genomic era. Nature 2000, 405(6788):823-826. 

63. Assenov Y, Ramirez F, Schelhorn SE, Lengauer T, Albrecht M: Computing 
topological parameters of biological networks. Bioinformatics 2008, 
24(Suppl 2):282-284. 

64. King RW, Peters JM, Tugendreich S, Rolfe M, Hieter P, Kirschner MW: A 20S 
complex containing CDC27 and CDC16 catalyzes the mitosis-specific 
conjugation of ubiquitin to cyclin B. Cell 1995, 81:279-288. 

65. Aristarkhov A, Eytan E, Moghe A, Admon A, Hershko A, Ruderman JV: E2-C, 
a cyclin-selectiveubiquitin carrier protein required for the destruction of 
mitotic cyclins. Proc Natl Acad Sci 1996, 93:4294-4299. 



Fang et al. Proteome Science 201 1, 9:26 
http://www.proteomesci.eom/content/9/1/26 



Page 15 of 15 



66. Yu H, King RW, Peters JM, Ki rschner MW: Identification of a novel 
ubiquitin-conjugating enzyme involved in mitotic cyclin degradation. 
Cuff Biol 1 996, 6:455-466. 

67. Vernoux T, Autran D, Traas J: Developmental control of cell division 
patterns in the shoot apex. Plant Mol Biol 2000, 43:569-581. 

68. Holt LJ, Krutchinsky AN, Morgan DO: Positive feedback sharpens the 
anaphase switch. Natufe 2008, 454(Suppl 7202)353-7. 

69. Haizhong F, Qingguo C, Jian F, Jian Z, Xiaohui Y, Jianru Z: Functional 
characterization of the Arabidopsis eukaryotic translation initiation 
factor 5A-2 that Plays a crucial role in plant growth and development 
by regulating cell division, cell growth, and cell death. Plant Physiology 
2007, 144:1531-1545. 



doi:1 0.1 186/1477-5956-9-26 

Cite this article as: Fang et al:. Comparative proteomics analysis of 
proteins expressed in the 1-1 and 1-2 internodes of strawberry stolons. 

Prnteome Science 201 1 9:26. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 

Submit your manuscript at (^\ RinMfid r pnt ^\ 

www.biomedcentral.com/submit Blomea ^ enirai 



